SPhT: T04/128 



Modeles a 1 et 2 matrices et fonction generatrice de cylindres a bords 

bicolores. 

N. Orantin 1 
Service de Physique theorique, CEA/Saclay 
Orme des Merisiers F-91191 Gif-sur-Yvette Cedex, France 



Abstract: 

Je vais traiter dans ce rapport des modeles a une et deux matrices comme outil pour 
etudier les theories conformes avec bords. J'introduirai dans un premier temps certains 
resultats bien connus que le lecteur pourra trouver en grande partie dans les ouvrages de 
references que sont |3T| . |H0] et [Hj pour ce qui est des matrices aleatoires et j20j pour la 
geometrie algebrique. Pour ce faire, apres une rapide presentation historique des matrices 
aleatoires et de leurs applications, je m'interesserai tour a tour aux modeles a une et deux 
matrices, presentant leur formalisme, expression en termes de diagrammes de Feynmann et 
surfaces aleatoires ainsi que leur lien avec les theories conforme. Je conclurai sur le modele 
a deux matrices en expliquant la methode des equations de boucles. Je m'attacherai dans 
un second temps a presenter un travail de recherche personnel consistant en le calcul de la 
fonction de correlation generatrice des surfaces aleatoires ayant la forme d'un cylindre aux 
bords bicolores. 

1 Introduction de la notion de matrices aleatoires et 
de leurs applications. 

Parmi les nombreuses applications des matrices aleatoires en physique que j'aborderai 
plus bas, je me suis particulierement interesse a leur interpretation en termes de surfaces 
aleatoires. En effet, les integrates de matrices peuvent etre vues comme les fonctions 
generatrices de surfaces discretisees avec ou sans bords. Les operateurs vivant sur ces 
surfaces sont decrit par les theories conformes dans la limite oil les surfaces discretisees 
deviennent continues. Les exposants critiques de tels operateurs sont bien connus lorsque 
ceux ci vivent a l'interieur de la surface ([22]). Les operateurs agissant sur les bords sont 
quant a eux bien moins connus et les modeles de matrices semblent etre un bon outil pour 
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determiner leurs exposants critiques. Les fonctions generatrices des surfaces avec deux et 
quatre operateurs sur un meme bord ayant deja ete determinees ([H], [TB] ) , ces notes ont 
pour but de presenter le calcul de la fonction generatrice des surfaces avec deux bords 
portant chacun deux operateurs. 

1.1 Historique. 

Le terme aleatoire associe aux matrices dont je vais introduire l'etude dans ce rapport sem- 
blerait induire que ces dernieres sont apparues comme un outil pour la physique statistique. 
Cependant, il en a ete tout autrement. 

En effet, elles ont vu le jour dans le cadre de la physique atomique et plus precisement 
dans le cadre de l'etude de noyaux lourds par Wigner en 1951 ([HI])- Le spectre en energie 
de ces derniers etant trop dense pour pouvoir tenir compte separement de chaque niveau, 
Wigner a considere le probleme d'un point de vue statistique. II a observe que la dis- 
tribution moyenne des niveaux d'energie, des fonctions de correlation ainsi que d'autres 
parametres tels que l'espacement moyen entre les niveaux d'energie, correspondaient au 
spectre d'une matrice aleatoire suivant une loi de probabilite gaussienne prise dans un 
ensemble dependant du probleme considere (je reviendrai plus loin sur la correspondance 
entre ensemble de matrices et caracteristiques du probleme considere). La surprenante 
fidelite des resultats experimentaux avec cette description en termes de matrices imposa 
pour la premiere fois les matrices aleatoires comme un puissant outil mathematique de 
description de systemes physiques complexes. 

Des lors, les matrices aleatoires ont connu un succes sans cesses grandissant en appa- 
raissant dans de nombreux domaines de la physique pourtant eloignes les uns des autres. 
Avant d'aller plus avant dans une presentation technique de cet outil, je vais brievement in- 
troduire certains de ses principaux domaines d'application ainsi que l'une de ses principales 
caracteristiques: le principe d'universalite. 

1.2 Le principe d'universalite. 

On peut considerer la notion de matrice aleatoire comme une generalisation des nombres 
aleatoires. On est alors en droit de se demander s'il existe un equivalent au theoreme 
centrale-limite dans le cas des matrices. Ce comportement particulier aux grands nombres 
peut etre observe lorsque l'on fait tendre la taille des matrices vers l'infini. Le spectre 
obtenu dans cette limite ne depend plus de la loi de probabilite initiale mais simplement 
des symetries de l'ensemble de matrices considere. 

Ce comportement n'est, pour l'instant, que conjecture. Cependant, les verifications 
experimentales ainsi que des demonstrations dans des cas precis permettent d'avoir suff- 
isamment confiance en cette propriete pour en faire l'une des principales composantes du 
domaine. Elle permet en effet de connaitre le comportement de n'importe quel ensemble 
de matrices, a grand nombre de variables, en utilisant une loi de probabilite gaussienne. 

II est alors important de classer les differents ensembles de matrices selon leurs pro- 
prietes de symetries. On rencontre trois categories de tels ensembles selon le point de 
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vue que Ton adopte, mais je n'en presenterai ici que deux, la troisieme correspondant aux 
matrices de transfert. D'une part, on peut considerer le hamiltonien comme une matrice 
aleatoire suivant une loi de probabilite gaussienne. On devra alors considerer les ensembles 
suivants: 

• Lorsqu'il n'y a aucune symetrie particuliere, on considere l'ensemble des matrices 
hermitiennes, appele Gaussian Unitary Ensemble (GUE), car il est invariant par les 
transformations unitaires; 

• Lorsque Ton ajoute une invariance du systeme par renversement du temps, il con- 
vient alors de considerer l'ensemble des matrices reelles symetriques, appele Gaussian 
Orthogonal Ensemble (GOE); 

• Si, de plus, le systeme a un spin total demi-entier et que la symetrie de rotation est 
brisee, on utilisera l'ensemble des matrices quaternioniques self-duales reelles, appele 
Gaussian Symplectic Ensemble (GSE). 

D'autre part, on peut considerer l'operateur d'evolution comme une matrice aleatoire 
suivant une loi de probabilite constante, ce qui fait apparaitre les ensembles suivants: 

• Sans symetrie particuliere, on aura l'ensemble des matrices unitaire ou Circular Uni- 
tary Ensemble (CUE); 

• Avec renversement du temps, on prendra les matrices orthogonales, soit le Circular 
Orthogonal Ensemble (COE); 

• Enfin, on peut considerer l'ensemble des matrices symplectiques, le Circular Sym- 
plectic Ensemble (CSE). 

1.3 Quelques applications. 

L'idee la plus naturelle d'utilisation des matrices aleatoires consiste simplement a considerer 
un systeme desordonne (|25j) ou le hamiltonien est lui-meme une matrice aleatoire. Selon 
les symetries precedemment enumerees, on aura le hamiltonien dans les ensembles GUE, 
GOE ou GSE. 

Une seconde application intervient lors de l'etude d'un probleme tout aussi naturel: le 
chaos quantique Si l'etude des systemes chaotiques est bien connue en mecanique 

classique, sa generalisation en mecanique quantique est tres problematique. La encore, 
les matrices aleatoires apparaissent indirectement. En effet, le hamiltonien du systeme 
considere n'est pas a proprement parler aleatoire. Cependant, les observations ont montre 
que, si le probleme classique associe est chaotique, son spectre est le meme que celui d'une 
matrice aleatoire appartenant aux ensembles GUE, GOE ou GSE. On peut egalement 
considerer l'operateur d'evolution qui, lui, appartiendra a l'un des ensembles CUE, COE 
ou CSE. Notons que cette relation entre chaos quantique et matrices aleatoires n'est pas 
demontree mais simplement observee experimentalement et conjecturee theoriquement. 
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Les matrices aleatoires apparaissent egalement en chromodynamique quantique (QCD). 
Cette theorie suppose l'existence de particules fondamentales constitutives des hadron, les 
quarks et les gluons, qui caracterisent l'echange de la charge de couleur (de dimension N=3) 
entre les quarks. Cet echange se traduit par une matrice 3x3. De maniere analogue a la 
mecanique quantique oil Ton doit integrer sur tous les etats possibles, il est ici necessaire 
de considerer un ensemble de matrices pondere par une loi de probabilite. L'utilisation 
du modele de matrices aleatoires a permis d'obtenir des predictions non perturbatives. 
Pour ce faire,il a ete propose par 't Hooft ([23|])de faire un developpement en puissance 
de 1/N quand le nombre de couleur N est artificiellement amene a tendre vers +00. Ce 
developpement en puissances de 1/N est en fait un developpement topologique sur lequel 
nous reviendrons plus loin. On peut egalement profiter du principe d'universalite en faisant 
tendre le nombre de degres de libertes du systeme vers l'infmi ([HI]). Enfin, cet outil 
mathematique permet de faire un lien entre la QCD et la theorie des cordes ([TT]. [TU] . 
23 , [24j) , sujet que nous allons aborder dans le prochain paragraphe. 

Les matrices aleatoires et la theorie des cordes se sont rencontrees a plusieurs reprises 
([TT]. [T2], [IH]). Je ne vais evoquer ici qu'un aspect de cette relation. Comme nous le 
verrons plus loin, les matrices aleatoires permettent, grace au developpement topologique 
mis en place pour la QCD, de decrire et de sommer des surfaces aleatoires. Or, les cordes 
evoluant dans le temps, forment justement de telles surfaces (pj, [H]). II est done naturel 
d'utiliser un modele de matrices pour traiter un tel probleme. Bien que les modeles de 
matrices permettent d'avoir une theorie de la sommation des surfaces, ils sont limites a 
l'etude de la theorie des cordes en dimension D < 1, ce qui empeche d'atteindre le cas 
physique D = 10, 26. Cependant, les modeles de matrices restent interessants que ce soit 
comme toy-model pour D < 1 ou grace aux autres liens qu'ils entretiennent avec la theorie 
des cordes et les recentes tentatives d' unifications comme la theorie M ([12], jl], EHj)- 

Au cours de cet expose, je reviendrai plus longuement sur deux autres importantes 
applications: l'etude des surfaces aleatoires et celle des theories conformes. 

2 Modele a une matrice. 
2.1 Introduction des notations. 

En premier lieu, introduisons les notations utilisees tout au long de ces notes ainsi que 
certaines considerations generales. 

Soit un ensemble E de matrices de taille NxN sur lequel on defini une loi de probabilite 
P(M) comme un poids de Boltzmann: 

P(M) = ^e-^ TrV ^dM (1) 
oil le potentiel V(M), la mesure dM (dans le cas ou E est l'ensemble des matrices 
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hermitiennes) et la fonction de partition Z sont dermis par: 



dM = Ui=i dM a Ui<j dReMij dlmM^ 
Z = J E dMe-^rV(M) =e -^F 

V(x) = V + 



(2) 



Note : le cas gaussien evoque precedemment correspond au potentiel V(x) = |i 2 . 



Nous avons introduit ici l'energie libre F qui est une fonction fondamentale dans l'etude 
des matrices aleatoires. On montrera ainsi facilement que : 



oil Ton a note (l\x l ) = (n^i)c+Ei(^i)(nj^£j}c+ + Ei^ 

Notons que les derivees secondes de l'energie libre par rapport aux ti correspondent a 
des fonctions de correlation a deux points et sont en fait des fonctions universelles. De 
maniere plus generale, toute derivee seconde de F est une telle fonction universelle ([E], 
[2]) et ce sont ces dernieres qui sont porteuses du principe d'universalite. 

2.2 Expression en termes de diagrammes de Feynman et de sur- 
faces aleatoires. 

Comme il Fa ete dit plus tot, l'une des principales applications des matrices aleatoires se 
rapporte a leur interpretation en termes de diagrammes de Feynman (jHj, PD]> [7j). 

Considerons l'ensemble E des matrices hermitiennes M. On aimerait calculer la fonction 
de partition Z definie plus haut. Afin d'obtenir une interpretation des modeles de matrices 
en gravite quantique, on va proceder par un developpement perturbatif autour d'un poten- 
tiel gaussien. En effet, on sait effectuer ce calcul dans le cas d'un potentiel gaussien et Ton 
pourra ainsi obtenir le cas general comme perturbation de l'integrale gaussienne. Notons 
l'analogie avec la physique des particules ou le potentiel gaussien correspond a des particules 
libres sans interaction, les integrates perturbees etant obtenues comme un developpement 
en diagrammes de Feynman. Dans le cas des matrices aleatoires, ce developpement ana- 
logue en diagrammes s'identifiera a un developpement topologique en termes de surfaces. 

Defmissons done 8V par : 




(3) 



k\...k. 



dt kl ...dt k . 





(4) 
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En considerant la partie quadratique au voisinage de son minimum comme dominante 
et effectuant le developpement a 5V petit, on peut ecrire le developpement : 

Z = [ dMe~^ TrM \l + ^Tr5V(M) + ^[TrSV(M)} 2 + ...) (5) 

Je 1 21 2 

Ce qui donne : 

Z = e- 2 *'/ E dMe-*™ 2 (e*™ l ') o 

= e -*> liMe-V**" (l + * (TrSV) + * {(TrW)\ + ...) (6) 

ou l'indice indique que l'on moyenne sur la distribution gaussienne / dMe~^~ TrM2 . 

Des lors, le theoreme de Wick permet de decomposer les moyennes de produits de vari- 
ables aleatoires en une somme sur tous les appariements du produit des valeurs moyennes 
par paires. 

Par exemple: 

(01020304) O = (0102> O (0304) O + {<j>lfa) (h<Pi} + (^Mo (0203) (7) 

C'est ce developpement qui est a la base de la representation diagrammatique de 
l'integrale. Pour representer chaque terme de ce developpement sous forme de diagramme, 
nous allons introduire des regies de Feynman expliquant comment les tracer et quel poids 
leur attribuer. 



2.2.1 Les regies de Feynman. 

Les regies de Feynman consistent a faire correspondre a chaque terme sous forme d'une 
valeur moyenne, une representation diagrammatique composee de vertex relies entre eux 
par des propagateurs. 

La regie fondamentale consiste a definir le propagateur correspondant a la valeur 
moyenne du produit d'une paire de variables gaussiennes. Le propagateur est donne, 
comme il est d'usage en theorie des champs, par l'inverse de la partie quadratique : 

j k im 2 

Notons que, contrairement a la theorie classique des champs, le propagateur n'est pas 
simplement une ligne non orientee, mais est compose de doubles lignes orientees. 

Les vertex, representant les termes de la forme ^^TVM*^, sont quant a eux issus 
de {TrSV) . Chaque terme de degre k sera represents par un diagramme a k lignes ex- 
ternes sous la forme de lignes doubles orientees de fagon a pouvoir les connecter a des 
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propagateurs. On a alors, pour le terme de degres k du potentiel 




2.2.2 Determination du developpement diagrammatique. 

Les regies de Feynman en main, nous avons les outils pour determiner le developpement 
diagrammatique de Z. On a, en effet: 

1 , Nt k 



r j^_»£ TrM 2 ^V, L n k r Tk 



e t J E dMe 2t 1 rm 



= e n^r(-^™ fe r (io) 



On peut developper chaque terme de la forme (jl k ^(-UTrM*)™'^ en une somme 
sur les diagrammes a n k vertex a k lignes externes dont le poids est: 

1 / T \ n * ( Nt 



oil Q est le facteur de symetrie du diagramme, n k le nombre de vertex a k pattes externes, 
n v = Y,k n k le nombre total de vertex et 1 le nombre de boucles du diagramme. 

Z s'obtient done comme la somme de tous les diagrammes formes de n k vertex a k 
pattes, n p propagateurs et comprenant 1 boucles auxquels on a associe le poids : 

}_ N n v -n p +l T n v -n Pt n P ^^_ tk y k ^ 
^ k 

2.2.3 Lien diagramme-surface. 

Les diagrammes definis plus tot peuvent etre interpretes de maniere geometrique. En effet, 
a tout vertex a k pattes, on peut faire correspondre un polygone a k cotes tel que chaque 
cote soit orthogonal a un propagateur. Par exemple, un vertex a trois pattes donne un 
triangle et un vertex a 4 lignes externes donne un carre. 



On obtient ainsi une relation de dualite mettant en correspondance un diagramme 
avec une surface discretisee. Le developpement precedent peut alors etre vu comme une 
sommation sur les surfaces discretisees. Ainsi, un diagramme ne comprenant que des vertex 
a trois pattes sera le dual d'une surface decoupee selon des triangles. 

Peut-on alors avoir une expression directe de Z comme duale de (JT2J)? Ceci n'est 
possible qu'en separant les diagrammes non-connexes des diagrammes connexes. Pour ce 
faire, introduisons une nouvelle fonction de partition. 

Soit Z = InZ . Le developpement diagrammatique de Z s'obtient en restreignant la 
sommation obtenue pour Z sur les diagrammes connexes: 

Z= Mn v -n p +l T n v -n vt -n P ^_ tk y k ^ 

diagrammes connexes k 

Z peut alors etre interprets comme la fonction de partition d'un ensemble de surfaces 
aleatoires tirees avec un poids de Boltzmann. Notons les interpretations des parametres 
dermis precedemment: tk represente la fugacite des polygones a k cotes; t<i joue le role 
d'energie de liaison entre les polygones; n p represente le nombre d'aretes; n v est le nombre 
de polygones const it ut if s de la surface; 1 est le nombre de sommets; x — n v — n p + I, 
exposant de N, represente la caracteristique d'Euler Poincare de la surface. 

2.2.4 Developpement topologique en 1/N 2 . 

Le developpement en puissances de N effectue plus haut peut s'interpreter de maniere 
geometrique. Considerons en premier lieu le developpement diagrammatique. % = n v — 
rip + l, l'exposant de N, est en fait une propriete de la surface sur laquelle peut etre dessine 
le diagramme considere pour ne pas avoir de lignes qui se croisent. Par exemple, les 
diagrammes correspondants a x = 2 peuvent etre traces sur un plan, x = sur un t° re 
etc... De maniere generale, si l'on note x = 2 — 2h, h est le genre d'une telle surface, c'est 
a dire, son nombre de trous. Ainsi, un diagramme de caracteristique x — 2 — 2h ne peut 
etre dessine que sur une surface avec au moins h trous. 

De meme, dans la description en termes de surfaces discretisees, le developpement de 
Z peut etre vu comme un developpement topologique. Ainsi: 

+oo 

Z = J2N 2 ~ 2h Z h (15) 

h=0 

oil Zh est la fonction de partition restreinte aux surfaces polygonales de genre h. On 
a ainsi un developpement topologique selon le genre des surfaces, qui donne acces aux 
fonctions Zh de genre h fixe. Notons que la limite N grand selectionne les surfaces de genre 



S 



0. On l'appelle limite planaire ou limite spherique. 



Z = 




+ N 




+ N 




0(N~ 4 ) (16) 



2.2.5 Interpretation de (TrM n ). 

Remarquons que le calcul de Z ne fait intervenir que des surfaces fermees duales de dia- 
grammes eux aussi fermes. Comment tenir compte des surfaces ouvertes, i.e. avec bords? 

Dans ce paragraphe, nous allons introduire de telles surfaces en nous interessant aux 
fonctions generatrices : 



T n = (TrM n ) = - f dM (TrM n ) e ~ NTrV ^ 
Z J 



(17) 



Du fait de la presence du facteur TrM n , chaque surface intervenant dans le calcul de 
T n contient au moins un polygone a n cotes. Le reste du diagramme est done une surface 
a laquelle on a enleve un tel polygone. C'est done une surface avec un bord de longueur 
n, et Ton a la relation : 



-T 

- 1 n. 



n 



E 



surface de bord n 



N X T n k -n p f - 



n 

k 



-t k ) 



(18) 



Que devient le developpement topologique defini pour les surfaces sans bords? II est 
toujours valable en tenant compte de la modification dans la definition de la caracteristique 
d'Euler Poincare x d'une surface avec un bord : x — 1 — 2/i. 

On peut alors ecrire le developpement topologique : 



T n = Y.N l ~ 2h T 1 



n[h] 



N 



+N 




+ ... 



(19) 



oil T n [h] est la fonction de partition partielle des surfaces de genre h ayant un bord de 
longueur n. 

II sera tres utile pour la suite de definir la resolvante qui jouera le role de fonction de 
partition grand canonique : 



1 +^ T n _1_+^ iTtM 



n=0 



+ 1 



N \ z- 



M 



(20) 



On peut generaliser la notion de surface avec bord en considerant les surfaces avec 
plusieurs bords. Pour ce faire, introduisons la fonction generatrice des surfaces ayant m 
bords de longueurs l±, l m definie par : 

T h ,..., lm = (TrM h ...TrM 1 ™) (21) 



9 



L 'equation analogue de (fTHjl s'ecrit alors: 

1 T h _ lm = £ ±N*T^ n n- 2 np m-tKT k (22) 



"' m surfaces de bords li, ...l m k 

ou l'on redefini \ = 2 — 2h — m. 

2.3 Matrices aleatoires et theorie conforme. 

Dans la partie precedente, nous avons presente le lien entre matrices aleatoires et surfaces 
discretisees. Nous allons maintenant nous interesser au passage a la limite des surfaces 
continues et ainsi aborder le lien avec les exposants critiques et la theorie conforme. 

2.3.1 Limite continue. 

Pour obtenir une surface continue, il faut faire tendre le nombre de pieces formant la surface 
vers l'infmi tout en rendant leur taille infmitesimale. Ces comportements singuliers doivent 
cependant repondre a certaines contraintes, telles que garder l'aire de la surface ( i.e. le 
nombre de pieces multiplie par leur surface ) ainsi que la longueur des bords finies. 

Une telle limite est obtenue lorsque l'un des parametres du potentiel ( respect ivement 
plusieurs parametres ) s'approche d'une valeur critique ( resp. d'un multiplet de valeurs 
critiques ). On est alors a un point critique ( resp. multicritique ). 

Soit n k le nombre de polygones a k cotes. En considerant l'expression de Z dans ()14|). 
on peut ecrire : 

. . dlnZ , , 

(n k ) = t k —— (23) 
atk 

En un point critique t^ c ou (n^) — + oo, la fonction de Z est singuliere avec un exposant 
critique a : 

% Z re g -\- Z s i n g (24) 

oil Z reg est reguliere et Z sing oc (tk — tk c ) a - Le nombre moyen de polygones a k cotes se 
comporte alors : 

(n k ) ex (t k - t kc ) a - 1 (25) 

Les exposants a sont universels et decrits par la theorie des champs conforme. Leur 
determination represente Fun des objectifs principaux des matrices aleatoires appliquees 
aux theories conformes. 

2.3.2 Action d'Einstein-Polyakov. 

En theorie des cordes et gravitation quantique, la fonction de partition s'ecrit sous la forme 

am): 

Z= J2 dSe- £ ^ (26) 
surface S 
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oil Paction d'Einstein-Polyakov est definie par 



£ = AnGx + A aire + matiere (27) 

G etant la constante gravitationnelle et A la constante cosmologique du modele. 
Supposant que les polygones constitutifs de la surface aient des cotes de longueur e, par 
identification de (|2*7jl avec ((Bj) , on obtient : 

InN = —AnG . , 

A k e 2 = ln(-t k ) - \ln(t 2 ) ^ } 

Dans ces conditions, l'aire peut s'ecrire A = (n t ) e 2 ou n t = n k est le nombre total 
de polygones constituants la surface. Alors (nt) doit diverger en \. Ceci nous amene a 
definir une constante cosmologique renormalisee A^ telle que : 



2 



t k = t kc - A kR e «-i (29) 
et Ton a l'ordre de grandeur de l'aire moyenne : A oc aA k ^ . 

Note : En fait, l'exposant critique a introduit ici, n'est pas bien defini. En effet, un tel 
exposant, valable pour tous les genres, n'existe pas. On peut seulement definir un exposant 
ah pour chaque ordre du developpement topologique. Cependant, le passage a la limite 
continue est toujours possible en procedant ordre par ordre en ■h. 

2.3.3 Theorie conforme et modeles minimaux. 

II existe tres peu de modeles que Ton sache resoudre completement ([22j): ce sont les 
modeles minimaux (p,q) dont l'une des caracteristiques principales, la charge centrale c, 
s'ecrit : 

(p ~ q) 2 

c = 1 — 6 avec p et q premiers entre eux. (30) 

pq 

Les dimensions des champs primaires sont alors les poids conformes de la table de Kac : 

h rs = - ^ - - ^ avec < r <p et < s <q (31) 
Apq 

Pour pouvoir comparer les resultats obtenus par les methodes de matrices aleatoires 
et de theories conformes, il faut se ramener a des quantites calculables dans les deux cas. 
Ce seront les exposants critiques. Parmi ceux-ci, les plus couramment utilises sont les 
exposants lies a l'aire et a la topologie, c'est a dire le genre h. 

On peut, par exemple, considerer le comportement, pour A — > oo, de la fonction de 
partition a aire et genre fixes : 

Z h (A) ~^oo A^- 3 (32) 
Ceci nous permettra de definir la susceptibilite de corde 7 s4r : 

7 A = 2-(2- 7atP )(l-/i) (33) 
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Dans le cas des modeles minimaux, en genre h = 0, on peut montrer que Ton a 



Istr = -2 J P , 9 ' (34) 
p + q — \p — q\ 

Remarquons que, les modeles de matrice donnant 7 str = — p+ ^_ 1 , les resultats ne coin- 
cident que pour les modeles unitaires. 

On peut egalement s'interesser au comportement d'autres operateurs O rjS vivant sur la 
surface : 

(O r , s ) (A) ~ A 1 '^ (35) 

L'exposant A rjS peut etre determiner par les theories conformes si l'operateur considere 
vit a l'interieur de la surface. Si c'est un operateur de bord, la methode des matrices semble 
etre une methode efficace en l'absence de theorie generate.. 

2.3.4 Exemple: cas d'une surface a deux bords; le cylindre. 

Dans cette partie, nous allons etudier le cas des surfaces a deux bords, decoupees selon 
des carres, comme exemple d'application. Notons qu'une surface a deux bords n'est rien 
d'autre qu'un cylindre. 

Nous allons done etudier la fonction generatrice des surfaces a deux bords, l'un de 
perimetre L = le et l'autre K = ke, ou e represente la taille du cote d'un triangle : 

T lih = (TrM l TrM k ) (36) 



c 



En utilisant le potentiel quartique symetrique: 

V(M) = X -M 2 - |M 4 (37) 
Comme on l'a vu precedemment, le nombre moyen de carres de cote e est : 

Ainsi, l'aire moyenne d'un cylindre de bords de perimetres L — le et L' — I'e sera : 

A = (39) 

II faut done determiner les fonctions T2i^k{g)- Pour ce faire, considerons maintenant la 
fonction a deux points, generalisation de la resolvante u>(z) : 

/ i i \ OO OO rp 

1 } \ x — M y-M/ c toto xl+1 y k+1 
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Dans la litterature ([HI), on trouve l'expression de u(x,y) qui est en fait le noyau de 
Bergmann (cf le paragraphe qui lui est consacre): 



u(x,y) 



1 



2{x-yf 



1 - 



xy — a 
a(x)a(y) 



(41) 



ou Ton a considere, profitant de la symetrie du potentiel, que la densite p(x) a pour 
support le segment [—a, a] (c'est ce que Ton appelle le cas a une coupure). On a defini la 



fonction a(x) = x 2 — a 2 . 



Par etude dimensionnelle, on sait que T k j = Ck^a k+l ou C^i est un nombre independant 
de g. On a ainsi : 



A(L, K) = e 2 (l + k)g 



dlna 
dg 



(42) 



On peut par ailleurs ([H3) determiner a 2 = ^- (l — y/1 — I2g) . Ce qui nous donne : 



A(L, K) = e(L + K)- 



2 1-^1-12^ 

3 9 



2g (l - v^!2i 



(43) 



Le passage a la limite continue necessite de considerer e — > tout en gardant L, K et 
A finis. Pour cela, on doit se placer au voisinage d'un point critique g — > g c = avec 
9 — 9c ~ e 5 ''■ Soit alors la constante cosmologique renormalisee A = -h, telle que : 



9 



12 



1-At 



(44) 



On obtient alors : 

A(L, K) = 36(L + K)Lq + 0(e) (45) 

Notons que ce resultat correspond a celui un cylindre de hauteur L , ce qui nous 
renseigne sur la forme de la surface. 



3 Le modele a deux matrices. 
3.1 Introduction du modele. 

Ce modele est une extension du modele a une matrice dans le sens ou l'on considere 
maintenant deux matrices hermitiennes Mi et M2 de taille N x N, ainsi que l'integrale : 

Z = J dM 1 dM 2 e -XTr{V,{M,)+V 2 {M,)-cM,M 2 ) (4g) 

avec Vi et V2 deux potentiels polynomiaux definis de maniere analogue au V precedent. 

Tout d'abord, on peut remarquer que l'integrale est modifiee de fagon triviale par 
changement de variable affine sur les matrices. Ainsi, on peut se placer dans le cas ou les 
coefficients des termes quadratiques de V\ et V2 sont tous deux eg 
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Des lors, une interpretation sous forme de diagramme peut etre introduite comme dans 
le cas du modele a une matrice, en prenant soin toutefois de differencier les deux matrices. 
On peut associer deux couleurs differentes aux deux matrices, bleu pour Mi et rouge pour 
M%. Les diagrammes de Feynman seront alors legerement plus compliques. 

lis feront ainsi apparaitre trois types de propagateurs differents : Tun tout rouge cor- 
respondant au lien M 2 -M 2 , l'autre tout bleu pour le lien M\-M\ et enfin un propagateur 
mixte avec un cote bleu et un cote rouge pour representer le lien entre M\ et M 2 . Notons 
que le poids de ce troisieme propagateur est egal a celui d'un propagateur unicolore multi- 
plie par c. Par contre, les vertex resteront quant a eux unicolores, aucun terme de la forme 
M® M| avec le produit ab > 1 n'apparaissant dans Taction. 

3.2 Potentiel cubique et modele d'Ising. 

Considerons le cas ou les deux potentiels sont cubiques. C'est a dire, f4*6j) s'ecrit : 

Z = JdM, dM 2 e -iVTr(f Mf+f Ml+lMf+iMl-cMiM.) (4?) 

Les diagrammes de Feynman sont alors constitues de deux vertex a trois pattes, l'un 
bleu et l'autre rouge, et des trois propagateurs definis precedemment : 




(48) 

Le passage aux surfaces duales ne fait plus alors apparaitre de simples surfaces trian- 
gulees, mais des surfaces composees de deux types de triangles : les uns rouges et les autres 
bleus. On peut representer ces differentes couleurs en attribuant par exemple un spin +1 
aux triangle rouges et un spin -1 aux triangles bleus. On identifie alors le modele a deux 
matrices a l'etude de surfaces triangulees sur lesquels sont disposes des spins + et -. 




(49) 



On note ainsi l'analogie avec le modele d'Ising. II reste a definir le lien entre les 
parametres de ces deux modeles (|2I|,[H!, 0) L 28j ) . 

Pour cela, il suffit d'identifier le poids associe a un diagramme dans chacun des deux 
modeles. 
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Pour le modele de matrices: un diagramme a n a triangles de spin a, n + _ propagateurs 

l p = 2' 



liant deux triangles de spins differents et n p = | {n + + n_) propagateurs en tout aura 



pour poids : 



N^gT 92~ c n+ ~(l - c 2 y n * (50) 



• Pour un modele d'Ising avec une constante cosmologique g, une constante gravitation- 
nelle N, un nombre total de triangle n, un champ magnetique h et une temperature 
J : 

N x ^g n e- J ^ a ' (7j e- h ^ a ' (51) 

En remarquant que 2n+_ = n p — J^^i^j, on peut identifier: 

, 3 



g = vgm 

e ~2h _ 31 ' (52) 

E 2J = c 



3.3 Methode des equations de boucles et applications. 

Dans cette partie, nous allons nous interesser a une methode de calcul permettant en 
partie de resoudre le modele a 2 matrices : la methode des boucles (JU21, [IH], |H], [I]), 
qui n'est autre que la methode de Schwinger- Dyson dans le cadre des matrices aleatoires. 
Dans un premier temps, nous aborderons cette methode d'un point de vue theorique avant 
de l'appliquer au calcul de la fonction ^ Tr x } Mi y \ l2 Tr x ,} Mi y il M2 ) Q ue j'ai effectue au 
cours de ce stage. 



3.4 Methode des equations de boucles. 

Considerons deux potentiels polynomiaux V\ et V2 definis par : 

Vl {x) = V p±L x k+l et V 2 (y) = V ^-y^ 1 (53) 

La methode des equations de boucles consiste simplement a effectuer un changement 
de variable infinitesimal de la forme : 

M i ^M i = M i + ef(M 1 ,M 2 ) pour i = 1,2 (54) 

dans les integrates de matrice. Notons que f(Mi,M 2 ) doit etre hermitienne pour 

que Mi le soit aussi. On va simplement ecrire l'invariance a l'ordre 1 en e de la fonction 
de partition Z ( on considere ici le cas % = 1, l'autre cas etant totalement analogue) : 

Z = J dM 1 dM 2 e- NTr{Vl{Ml)+V2( - M2) ~ MlM2 ' ) 

= J dMi (IM2 e -^i(Mi)+V 2 (M a )-M 1 M 2 ) 
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= J dM x dM 2 (1 + e J(M 1 ,M 2 ))e~ NTr ( Vl ^ Il)+v ^ M2) - A ~ IlM2 ) (1 - eK (M u M 2 )) + <3(e 2 ) (55) 

oil J (Mi, M 2 ) et K(Mi, M 2 ) correspondent a l'ordre 1 en e respectivement du jacobien 
du changement de variable et de la variation de Faction. C'est a dire : 

det Ht = 1 + eJ ( M l' M 2)+0(6 2 ) 

NTr (Vi(Mi) + V 2 (M 2 ) - MiM 2 ) = NTr (Vi(Mi) + V 2 (M 2 ) - M X M 2 ) 

+eK{M 1} M 2 ) + 0(e 2 ) 

(56) 

On obtient ainsi l'equation de boucle generique : 

(J(M 1 ,M 2 )) = (K(M 1 ,M 2 )) (57) 

3.4.1 Determination de J(M 1 ,M 2 ) et de K(M U M 2 ). 

• La variation de Paction par ce changement de variable donne aisement l'expression : 

K(M U M 2 ) = N {Tr [V{{M 1 )f{M ll M 2 ) - M 2 f(M h M 2 )]} (58) 

• La determination de J (Mi, M 2 ) est quant a elle moins evidente, mais peut heureuse- 
ment etre resumee par deux regies simples, connues sous le noms de regies "split" 
et "merge" (|14J). Celles-ci correspondent en fait aux deux types de changement de 
variable effectues en pratique. 

Regie Split. 

Le premier type de changement de variable usuel correspond a /(Mi, M 2 ) = A x J~ M B 
oil A et B sont des fonctions de Mi et M 2 . Alors la correction issue du Jacobien est : 

J(Mi,M 2 ) = Tr(A 1 —)Tr( ^-b] 

+ contributions venant de A(M l ) et B(M X ). (59) 

Ainsi, chaque fois que Ton rencontre un terme x \ lx en dehors d'une trace, on "coupe" 
l'expression en deux traces en introduisant un facteur x _ M dans chaque trace. 

Regie merge. 

On rencontre egalement des changements de variable du type 
/(Mi, M 2 ) = ATr (j^B). On a alors : 

J(M l5 M 2 ) = Tr(A 1 —B -—] 

K ' ' \ x-M x x-Mj 
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+ contributions venant de A(M l ) et B(M 1 ). (60) 

C'est a dire que, chaque fois que l'on rencontre un terme x J Ml dans une trace, on 
regroupe toute l'expression a l'interieur d'une meme trace en remplagant "Tr" par 
un duplicata du facteur X 2 M ■ 

3.4.2 Quelques definitions. 

Les equations de boucles faisant intervenir de nombreuses expressions fastidieuses a ecrire, 
il est necessaire d'introduire quelques notations pour que ce rapport soit lisible. Nous 
considererons done les fonctions suivantes : 

W(x) = ^-(tt ?— - \ et W(y) = ^-(Tr l — \ (61) 

v ' N\ x-Mxl yy > N\ y-M 2 / 

m . 1 / m 1 1 \ „, v 1 / m V{(x) - W(Mi) V4(y) - V4(M : 

W(x, y) = — ( Tr — — > et P(x, y) = — ( Tr ' 2yy> : 

v 7V\ X -Miy-M 2 / y ' N\ x - M x y - M 2 



iv\ x^ ~y~^~M 2 / x~^M l J^M 2> 

(63) 

(64) 
(65) 



1 / 1 1 1 

^ " * ) = ^ ( Tr ^ ^M 2 Tr 

£/(x v- x') = — ( Tr - Yi^JL^MlL T r _ 



L(x, y; y') = / Tr l — — Tr -V{{M{)\ (66) 

v ,y,yj N 2\ X -M iy -M 2 y'-M 2 1V V V ' 

v > )= 1 / Tr ^i / (^)-K(Mi)y 2 / (y)-y 2 / (M 2 ) 1 1 

* ' N 2 \ x - Mi y - M 2 x' - Mi y> - M 2 , 



, , _ 1 / 1_ _1_ 1 1 

H(x,y,x,y) - — ^ Tr - _ — - _ — - _ — - _ - 

1/1 1 1 1 

+ 2iV\ r x - Mi y' - M 2 x' - M x y - M 2 / ^ ' 

F{x,y,x,y) - — ^ Tr — — ——————— 

1 /_ V!(x) - VUM-i ) 1 1 1 



1 (r lY Vii.x)-V[{M 1 ) 1 
2iV \ x - Mi y> - M 2 x' - M x y - M 2 



(69) 



^ y; y,) = ^ ( Tr x^mT ^m 2 Tr x^mT T^m, ) (70) 

U( X y x > y >) = J_ / Tr ^)-^(Mi) 1 _J 1_\ 

1 ' y ' ' yj iV 2 \ x-Mi y-M 2 x'-M iy '-M 2 ) {<} 

Nous aurons egalement besoin de deux fonctions fondamentales : 

X(y) = V 2 \y)-W(y) et Y{x) = V((x) - W{x) (72) 
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3.4.3 L'equation de boucle maitresse. 

Parmis les nombreuses equations de boucles que l'on peut former, l'une a une signification 
particuliere et on l'appelle "master loop equation" (jS2j, IE])- 
Considerons le changement de variable : 

M 2 -> M 2 + e l -— (73) 

x-M 1 v 7 

L'equation de boucle associee est alors simplement : 

xW(x) - 1 = i ^Tr-^-y 2 '(M 2 )^ (74) 
Le changement de variable (hermitien bien stir) : 

Mi — > Mi + t( 1 - ^> - ^ + - ^> 1 „ ) (75) 
^ar-Mi y-M 2 y - M 2 x-M x ) K ' 

donne, en utilisant (JHJ), l'equation de boucle : 

(y - Y(x))U(x, y) = V 2 '(y)W(x) - P(x, y) - xW(x) + 1 - -^U(x, y; x) (76) 

U(x, y) etant un polynome en y, il n'a pas de singularity pour y fini, et l'equation (JZHJ) 
prise en y = Y(x) s'ecrit : 



V = Y(x) 

{V 2 \y) - x){V{{x) -y)- P(x, y) + 1 = &U(x, y; x) 



(77) 



En notant le polynome de degres {d\ + 1) en x et (d 2 + 1) en y : 

E(x,y) = (V 2 \y) - x){V{{x) - y) - P(x,y) + 1 (78) 
on peut reecrire l'equation precedente en : 

E(x,Y(x)) = j^U(x,Y(x);x) (79) 

qui est appelee equation de boucle maitresse. 

Cette equation est fondamentale car elle permet de determiner la fonction Y(x) et ainsi, 
donne acces a toutes les fonctions generatrices de surfaces aleatoires. 



Remarques: 

• En echangeant les indices 1 •*->• 2 dans tout le paragraphe precedent, on aurait obtenu 
une equation analogue a (|TH|) : 

E(X(y),y) = ^U(X(y),y;y) (80) 



18 



A l'ordre en e, les equations (|79*|l et (jSUj) nous donnent : 

( E(x,Y(x)) = . , 

\E(X(y),y) = ^ 

Les zeros de E(x,y) respect ivement en x et y seront alors notes X (y),Xi(y), X dl (y) 
et Y (x),Y 1 (x) i ...,Y d2 (x). D'ou : 



di d2 

E(x,y) = -g dl +i Y[(x - X t (y)) = -g d2 +i ]J(y - Y t (x)) 

i=0 i=0 



52) 



E(x,y) = defini une courbe algebrique. On voit ainsi apparaitre, pour la premiere 
fois le lien entre modeles de matrices et geometrie algebrique. 



3.5 Interpretation des traces et operateurs de bord. 

Comme nous l'avons vu plus haut, dans le modele a une matrice, les surfaces a k bords 

sont representees par les termes de la forme (^(Tr x } M ^j \ . 

Nous allons maintenant nous interesser a la generalisation de cette interpretation dans 
le modele a deux matrices. En effet, ce cas est bien plus complexe : du fait de la presence de 
deux matrices, on peut observer deux types de traces; certaines dites mixtes font apparaitre 
a la fois les matrices Mi et M 2 a l'interieur d'une meme trace, les autres etant dites non 
mixtes. 

II existe cependant une procedure simple pour interpreter les termes de la forme : 



Tr- 



x 



Mi 



n+ 



Tr- 



y 



Mo 



n 

k=l 



Tr 



x - Mty - M 2 



Les exposants introduits ici caracterisent les surfaces intervenant dans le developpement 
topologique. Ainsi, on pourrait montrer que : n + est le nombre de bords entierement 
constitues de segments +; n_ est le nombre de bords entierement constitues de segments 
-; n t = n + + n_ + J2k n k est le nombre total de bords;rifc est le nombre de bords constitues 
de la succession de k groupes de segments + contigus separes par k groupes de segments - 
contigus. Par exemple pour n + = 0, n_ = et n 2 = 1: 



Bord constitue de - 



Bord constitue de + 



Bord constitue de + 




(84) 



Bord constitue de 
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Ainsi, le nombre de traces indique le nombre de bords tandis qu'a l'interieur de chaque 
trace, le nombre de facteur x \ {i y \ l2 indique la moitie du nombre d'operateurs de bords, 
i.e. la moitie du nombre de fois ou le bord change de couleur. 

3.6 Derniers outils necessaires au calcul: le noyau de Bergmann 
et les operateurs d'insertion de boucle. 

3.6.1 Les operateurs d'insertion de boucles. 

On peut se demander comment passer d'une fonction a une boucle composee de +, W(x), 
a une fonction a deux boucles +, Vt(x]x') = (Tr x _ 1 Mi Tr x ,^ Mi ) ■ De maniere generale, 
on aimerait trouver un operateur qui insere une boucle (de type + ou - ) dans la surface 
considered. 

De tels operateurs existent (pQ, [2]). Considerons les derivations par rapport a V\(x) et 
V2(y) definies formellement par : 

d f d _ ^ k d _ 



9Vi(x) ti xk+ld 9k dV 2 (y) tiV k+ld 9k 

On peut alors remarquer que Q,(x]x') = Jp^-M • On pourrait egalement verifier que 
la derivation par rapport a V\{x) appliquee sur un terme de la forme a pour effet 
d'augmenter n + d'une unite. De meme, la derivation per rapport a ^(y) augmente n_ 
de 1. Ainsi, ces operateurs ont pour resultat l'insertion d'un bord (ou boucle) compose 
respectivement uniquement de + ou uniquement de -. Pour cela on les appelle operateurs 
d'insertion de boucle. 

3.6.2 Le noyau de Bergmann. 

Considerons a nouveau la fonction a deux boucles Q(x; x'). Elle peut s'ecrire: 

n(x-x') = ^l = l - (86) 

{ ' } dV l (x') {x-x'Y dV 1 {x') [ } 

Pour etudier le second terme du membre de droite, plagons nous sur la surface de 
Riemann definie par notre courbe algebrique £ = {(x,y)\E(x,y) = 0}. Un point p de £ 
est alors un couple de complexes p = (X(p),y(p)) tels que E(X(p),y(p)) = 0. 

Notons qu'a chaque x (resp. y) correspondent d 2 + 1 (resp. d\ + 1 ) points de £ : 

X(p) = x <^> p = pk{x) pour k allant de a d 2 - 
y{o) — V Q — Qk{ x ) pour k allant de a d\. (87) 
Des lors, on peut definir la fonction sur la surface de Riemann : 

B(p, P ') dy(p) 



dX(p)dX(pf) dV x {X{p')) 



X(p)=cte 



Cette fonction a un seul pole double en p=p' avec residu nul sur la surface de Riemann 
c'est le noyau de Bergmann ([20], qui s'appelle fonction de Weierstrass en genre 1. 
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4 Application au cas d'un cylindre a bords bicolores; 
Calcul de ( Tr — \j — \y- Tr -, * -, * ) . 

\ x-M 1 y-M2 x'-Miy'-M 2 / c 

Dans cette partie, je vais presenter ce qui fut l'essentiel de mon travail de recherche au cours 
de ce stage, c'est a dire le calcul de W(x,y;x',y'). Celui-ci permettra par la suite d'avoir 
des resultats relatifs aux cylindres dont chaque bords est compose de deux segments de 
couleur differente. 




(89) 



On pourra ainsi passer a la limite conforme pour obtenir des informations sur les ex- 
posants critiques en presence d'operateurs de bords. 

Pour effectuer ce calcul, nous allons proceder par etapes: 

• Determination d'un systeme d'equations de boucles faisant intervenir W(x,y;x',y'); 

• Resolution de ce systeme en commenceant par les inconnues dont les proprietes poly- 
nomiales rendent le calcul plus simple; 

• Simplification de l'expression de W(x,y;x',y') trouvee. 



4.1 Les equations de boucles. 

Dans cette partie, nous allons determiner deux equations de boucles dont les seules fonc- 
tions inconnues sont W(x,y;x',y'), U(x,y;x',y') et P(x,y;x',y'), les autres moments inter- 
venant ayant deja ete calcules dans la litterature. 

Notons que ces equations ne peuvent pas s'ecrire directement en appliquant la methode 
des boucles presentee plus haut. En effet, si usuellement seuls les termes dominants ( de 
plus haut degre en N ) s'annulent grace a une equation de boucle d'ordre inferieure, ici des 
termes sous-dominants disparaissent egalement car ils interviennent dans cette equation 
d'ordre inferieur. Ainsi, je preciserais pour chaque equation le changement de variable 
principal et celui permettant d'eliminer les termes dominants et certains sous-dominants. 

La premiere equation est obtenue en considerant le changement de variable: 



+ 



x - M x y - M 2 y- M 2 x-M 1 

i i 

x-Mi y-M 2 



L'equation de boucle associee a dM x = } M \ 4 + h.c, nous permet d'ecrire a l'ordre 



sous dominant : 

(Y(x) - y)W(x,y;x',y') - U(x,y';x',y') = W(x' ,y'; x)(W(x,y) 
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(91) 



6M> 



D' autre part, les changements de variables 
1 



V{{x)-V{{Mi) 1 1 V{{x)-V{{Mi) 

a: -Mi y-M 2 y - M 2 x - Mi 



Tr 



x' - M 1 y' - M 2 



(92) 



et <9M 2 = KMzgW i + fr. c donnent : 

* x— Mi y—M2 



(X(y) - x)U(x, y; x', y') - P(x, y; x', y') = L(x', y'; y) + U(x, y) - V{{x) W(x\ y'; y) 



(93) 



+ ; [F(x, y', x', y') - F(x, y, x', y')\ 

y-y 



4.2 Resolution des equations. 

Nous allons maintenant resoudre ces equations en utilisant le fait que P(x,y;x',y') est un 
polynome en x et y. 

Tout d'abord, une autre equation de boucle donnant : 



L(x\ y; y) = yW(x', y ; y) + H(x , y, x , y') 

,y)-Vl(x)+y 
[TE]' [T7j). on peut simplifier (|9*5|) en : 



et les resultats U(x, y) - V{(x) + y = ^ et W(x, y) = 1 - { J^_ Y) 



(94) 

etant connus 



(X(y)-x)U(x ) y ) x' ) y')-P(x,y;x',y') = ^l$W(x' >y ';y) 



(95) 



x-X 

+H{x', y, x',y') + — ^ [F(x, y', x 1 , y 1 ) - F(x, y, x', y')] 



Les differents termes du membre de droite peuvent etre calcules a partir de resultats 
preexistants (|IE|) : 



B=^lw(x',y';y) 
x — X 



E(x,y)E(x',y') 



x- X)(y> -Y>){x' - X>) 
^dX k (y>) 1 dY(x') 1 



fi dV *(y) X '- X 'k dV 2 (y)y>-Y> 



(96) 



et 



C = H(x', y, x', y') + [F(x, y', x', y') - F(x, y, x', y')) 

y-y 
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1 



(y-y')(x'-X')(y'-Y') 



E(x', y)E x (x', y') - E x (x' , y)E(x', y') 



+ 
+ 
+ 
+ 



E(x', y')E y (x, y 1 ) - E y (x', y')E(x, y') 



(x>-X)(y-Y>) 



(x - X')(x - x') 
E(x',y)E(x',y') 
(x'-X')(y-Y')(y'-Y') 
X-X' 



Y' 



(x> -X)(y-Y>)(y> -Y>) 
1 



(x> - XY(x> - X'){y - y>) (x - x'){y - y')(x> - X) 
E(', y ')E(x',y') 



(y-y>)(x-X>)(x> -X')(y'-Y') 

E(x', v r 



X' y 



(x-x'){y' -Y') {x - X'){x' - X') 
E(x,y)E(x',y') 



(x - x')(y - y')(x' - X') 2 (y' - Y') 2 {x - x')(y - y') 2 (x - X){x' - X')(y> - Y') 
, E(x',y)E(x,y') 

+ (x - x')(y - y') 2 {x - X')(x> - X)(y - Y') K *'> 



Notons que U (x, y; x', y') est un polynome en x. Ses poles en y ne dependent done pas 
de x. Alors (X(y) — x)U(x, y; x', y') s'annule pour tout y = Y^x). D' autre part, P(x,y;x',y') 
est un polynome en x et y de degre d 2 — 1 en y. II est done entierement determine par sa 
valeur en d 2 points, par exemple les Y^x) pour i = l..d 2 - 

Remarquant que E ^'x v = E x (x,Yi(x)), la formule d'interpolation prise pour les 
y = Yi(x) donne : 



y=Yi 



p( , ^ P(x,Y l -x\y')(Y,-Y) E(x,y) 

k (V-Y^y-Y) E y (x, Yi ) 

Soit, sachant que X' y = et Y< = -f^$ 



(98) 



d 2 



-P(x,y;x',y') = £ 



+ 



E x (x,Yj) E(x iy )E(x' jy ')(Y t -Y) 
- E y (x, Y) (y - Yi)(y - Y)(x> - X>)(y> - Y>) 

8X' k 1 dY' 1 

k[ dV2 ( Y i) x '- x 'k~ 9V 2 (Y) y> - Y> 

'_ (Yj-Y)E(x,y) 

(x - x')(Y - y>)(x> - X>)(y> - Y')(y - Y)(y - Y)E y (x, V,) 

E(x', y')E x (x', Yi) - E x (x>, y')E(x', Y) E(x', y')E y (x, y') - E y (x>, y')E(x, y') 



+ ■ 



(Yi Y') 
(Y i -Y)E(x',Y i )E(x,y)E(x',y') 



(x X') 



(x - x')(x' - X')(Y t - Y>)(y> - Y')(y - Y t )(y - Y)E y (x, Y) 
Yi x-X' 1 



(Yi — Y')(y' — Y') ( x - x>)(x> - X>)(Yi - y>) (x - x>)(Y - y>) 
(Y i -Y)E(x,y)E(x,y')E(x',y') 



(Yi 


-y')(x- 


-X')(x'-X')(y'-Y')(y- 


Y i )(y-Y)E v (x,Y i ) 




1 






[(x 


-x')(y> 


-Y') (x - X')(x' - X 1 ) _ 





(Y i -Y)E(x,y)E(x',y'f 



(x - x')(Y - y')(x' - X') 2 (y' - Y') 2 (y -Y t )(y- Y)E y (x, Y t ) 

(Y -Y)E(x,y)E(x> ,y')E x (x,Yi) 

(x - x')(Yi - y') 2 (x' - X')(y' - Y')(y - Y)(y - Y)E y (x, Y) 
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{Y-Y)E(x,y)E(x',Y)E(x,y') 

(x - x'f{y - y'f{x - X')(Y t - Y')(y - Y 2 ){y - Y)E y (x, Y % ) 



(99) 



Nous pouvons alors exprimer W(x,y;x',y') en fonction de P(x,y;x',y') en eliminant 
U(x,y;x',y') dans et <|55^ : 



(Y -y)W{x,y;x',y' 



X-x 

E(x,y) fW{x',y';y) W{x',y';xY 



x-X 



y-Y 



x-X 
+ H(x',y,x',y') 



X-x 

+7 — — y'i x> > y') - F ^ y> x '> f ')] ( 10 °) 

\y-y)\X-x) 

Utilisant les expressions de H(x,yx',y') et F(x,y,x',y') etablies dans [THj ou Ton fera 
tendre x vers x' et y vers y' ainsi que W{x', y'\ x) = 9 ^y^) ^ e ^ W{x', y'; y) = on 
obtient finalement l'expression de W(x,y;x',y') : 



E(x,y)E(x',y') 



(Y-y)W(x,y; x',y') 



1 



{x-X){x'-X'){y'-Y') 
, v d 2 E x (x,Yi) (Yj—Y) 



^k=l X '-X' 



1 



9XL 



1 dXL 



E y (x,Yi) (y—Yi)(y—Y) 



^k=l 



x-X dV 2 (y) y-Y dVi(x) 
1 dY' 1 



y'-Y' 



dV 2 (Y)x'-X' k dV 2 {Y)y'-Y' 
dY 1 1 dY' 



c-XdV 2 {y) y-YdVx{x) 



E(x,y)E{x' iy ') 



1 



1 {x-x'){y-y'){x-X){x'-X'){y'-Y') 

E{x,y')E{x',y') 

(y-y')(x-X)(x-X'){x'-X'){y'-Y') 

E{x',y)E(x',y') 

(x-x'){x-X) (x'-X') {y-Y' ) (y'-Y') 

E{x\y)E(x,y') 

{x-x'){y-y')(x-X'){x'-X)(y-Y') 

E(x',y' 



+ 



1 



(x-x')(y-Y) 1 (y-y')(x-X) 

i K 



(x-x'){y'-Y') (x-X')(x'-X') 

1 Yj 

(y-y')(x'-X>) {y-Y'){y'-Y') 

1 i 1 



+ 



(y-y')(x-X) 1 (x-x')(y'-Y) 



1 (x-x'){y-y'){x-X){x'-X') 2 (y'-Y') 2 
E(x',y)E x (x',y')-E x (x',y)E(x',y') E(x' ,y')E y {x iy ' )-E v {x' ,y')E{x,y') 



{x-x'){y-y'){y-Y'){x-X){x'-X'){y'-Y') (x-x'){y-y'){x-X){x-X'){x'-X')(y'-Y f ) 



J_ T u 2 



x-X 



(Yj-Y)E{x,y) 



{x-x'){Y t -y'){x'-X'){y'-Y'){y-YA{y-Y)E y {x,Yi) 



E(x',y')E x (x',Y)-E x (x',y')E(x',Yi) , E{x' ,y')E y {x,y')-E y (x' ,y')E{x,y' ) 
{Y-Y') "T {x-X') 

{Y-Y)E(x',Y)E{x,y)E{x',y') 

{x-x'){x'-X'){Y-Y'){y'-Y'){y-YA(y-Y)E y (x,YA 

YL x-X' , 1 



(Yi-Y>)(y'-Y>) (x-x'){x'-X')(Y-y') + (x-x'){Y-y') 
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(Y-Y)E(x,y)E(x,y')E(x',y') 



(Y-y') (x-X>)(x>-X>) (y'-Y') (y-Yi)(y-Y)E y (x,Yi] 



X' 



(x-x')(y'-Y') (x-X')(x'-X') 



(Y l —Y)E(x,y)E(x l ,y 



(x-x')(Y l -y')(x'-X') 2 (y'-Y') 2 (y-Y i )(y-Y)E y (x,Y) 

(Y-Y)E(x,y)E(x' ,y')E x (x,Yj) 

" (x-x')(Y l -y') 2 (x'-X')(y'-Y')(y-Y l )(y-Y)E y (x,Y) 

(Y l -Y)E(x,y)E(x',Y)E(x,y') 

(x-x') 2 (y-y') 2 (x-X')(Y l -Y')(y-Y l )(y-Y)E y (x,Y l ) 

(101) 



4.3 Simplification de l'expression de W(x,y;x',y'). 

L'expression obtenue plus haut permet un passage a la limite conforme. Cependant, du 
fait des nombreuses sommes sur les Y$, elle n'est pas tout a fait satisfaisante. On voudrait, 
en effet, trouver une expression qui montre que W(x, y; x', y') est une fonction definie sur 
£. Pour cela, il faut que cette derniere soit symetrique en les Y{. Nous allons done la 
simplifier. 

Tout d'abord, on peut noter que : 

^ Yj-Y 1 y-Y y'-Y \ 



\ (Yi - y')(y - YjE y (x, Y t ) y - y >^\(y- Y t )E y (x, Y t ) (y> - Yi)E y (x, Y { 

1 ( y-Y _ y'-Y ' 
y-y'\E(x,y) E(x,y>), 

Ceci nous permet de calculer trois sommes : 



(102) 



^2 (Y-Y)E(x,y)E(x',y') 2 



(x-x')(Y t -y')(x> -X') 2 (y> ' -Y') 2 (y-Y l )(y-Y)E v (xXi) 

E(x',y') 2 , (y'-Y)E(x' 7 y') 2 E(x,y) 

(x-x')(y-y')(x'-X') 2 (y'-Y') 2 ^ (x-x')(y-y')(x'-X') 2 (y'-Y') 2 (y-Y)E(x,y') 



(103) 



et 



(Yj-Y)E(x,y) 



(x-x')(Y i -y')(x'-X')(y'-Y')(y-Y i )(y-Y)E y (x,Y i ) 

1 (y'-Y)E(x,y) 



(x-x')(y-y')(x'-X')(y'-Y') (x-x')(y-y')(x'-X')(y'-Y')(y-Y)E(x,y') 



(104) 



v rf 2 (Y-Y)E(x,y)E(x,y')E(x',y') 

^»=1 (Y i -y')(x-X l )(x l -X l )(y'-Y')(y-Y i )(y-Y)E y (x,Y i ) 

E(x,y')E(x' ,y') (y'-Y)E(x,y)E(x' ,y') 



(y-y') (x-X>) (x'-X') (y'-Y') (y-y') (x-X>) (x'-X>) (y'-Y') (y-Y) 



(105) 
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Pour calculer les termes restants, on doit connaitre A = Y^=o (F^FTT^^Wfo y~) ■ P° ur 
ce faire, on va calculer A de deux manieres differentes. 
En ecrivant Yi — Y — Yi — Y' + Y' — Y, on obtient : 

A ! | E(x',y) * E(x', Yi ) 

1+ E(x,y) +{Y Y) h(Y*-Y')(y-Yi)E y (x,YA (1 ° 6) 

Alors que Yi — Y = Yi — y + y — Y donne : 

^-^-^ W-J-WT.) (107) 

Ceci nous permet d'etablir l'expression de A : 

y-Y 

y-Y' E(x,y) 

On calcule egalement : 



y-YE(x'y) , . 



^(y,-yO(y- W*-yO^M) y-y' 1 J 

c _y^ (y-y^yp _dB 

^ ) {Y t -Y'Y{y-Y l ){Y l -y')E y {x,Y,d SY' 1 J 

D f (y-y)£,«yp gg 

^(y-yO(y-yp(y-yO^,y) dx> * u 1 j 

* (y-y)£«yp &4 
^(y,-y0 2 (y-yp^(x,yp ay 1 J 

Le calcul des termes restants est assez different. En effet, ces derniers font intervenir 
des noyaux de Bergmann et ne dependent pas uniquement de E(x,y). 

Tout d'abord, exprimons EfcLi ( x >-x> k )dx av 2 { y ) ~ (y'-Y')dx Mk) 611 fonction des Y i et non 
plus des X k : 



et 



^ dX' k 1 ^ B(q,q' k ) 1 

h d Mv) * - K h dy{q)dy{q' k ) x(p>) - x(q' k ) 

= -I Res ^ B M 



(113) 



t^i y(t) - y' dy(q)dy(t) x(p') - x(t) 

Res jgft) 1 

9o.9^i y (*) - y' dy(q)dy(t) x{p) - x(t) 
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Ce qui peut s'ecrire : 



di 

E 



1 B(q' ,q) 



{dV 2 (y)(x'-X' k ) x'-X' dydy' (x> - X)(y - y'f 



D' autre part, on sait que : 



" Q (y> - YfidydxW 



B(p',q) 



(x-X)dV 2 (y) (x-X)dx'dy 
Et la meme methode donne les derniers termes : 



A dX' k 1 
h^x) (x'-X' k ) 



1 B(q' ,p) 



+ 



1 



x' — X' dxdy' (x' — x) 2 (Y — y') 
fl(pj,p) 



~"o (y' - Y{)dxdx{p' i ) 



(114) 



(115) 



(116) 



et 



di 

E 



9X1 



1 £(<7o, Pi) 



fe =i ^W) (x' - X' k ) x' - X' dy( Pl )dy' (x> - x)(y' - Ytf 



+ EtI7 



3 =o (v' - y;)dy{Pi)dy{p' j ) 
dY> B(p', Pi ) 



dV 2 (Y t ) dx'dy( Pt ) 



BY' 



B(p',p) 



(y-Y)dV 1 (x) (x-X)dx'dx 
Finalement, on obtient pour W(x,y;x',y') l'expression : 

(Y - y)W{x,y; x',y') = 

E(x,y)E(x',y') \ B(q',p) B(q\q) 



(x-X)(x'-X')(y'-Y') \ {y-Y){x'-X')dxdy' (x-X)(x> -X')dydy> 

7 + E 



(117) 

(118) 
(119) 



I ! \ I Btfn pj) _ B(p' Pi ) _,_ v d 2 B{p> pPi ) \ 



y-Y y-Yi J \(x'-X')dy'dx (y'-Y')dxdx' 1 ^j=0 (y'-Yfidxdyip'j) J 
1 ( B(p>,g) j_ B(p',,p) 



1 



y'-Y( \{x-X)dydx{p' i ) 1 [y-Y^xdxip'i) 



E{x,y)E{x' V) 



(ar-a;') (y-y') (x-X) (x'-X') (y'-Y>) 
(y-Y)E(x',y)E(x,y') 



1 



+ 



(x-x'){y'-Y) 1 (y-y')(^-*) 
_| E(x,y)E(x,y') 



{x-x') 2 (y-y') 2 {x-X')(x'-X)(y-Y')(y'-Y) ^ (x-x') 2 (y-y') 2 (x-X)(x-X')(y-Y) 
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(y'-Y)E(x,y) 



E(x',y')E v (x,y')-E v (x',y')E{x,y') 



(x-x')(y-y')(x-X)(x'-X')(y'-Y')(y-Y)E(x,y') 
(y'-Y)E(x iy )E(x',y') 



(x-X' 



{y-y'){x-X){x-X'){x'-X'){y'-Y'){y-Y) 



X' y 



(x-x')(y'-Y') (x-X>)(x'-X') 



(y'-Y)Y'E(x,y)E(x>,y>) 2 



(x-x')(y-y')(x-X)(x'-X')(y'-Y')3(y-Y)E(x,y') 

ylE{x' lV 'f 

(x-x')(y-y')(x-X)(x'-X')(y'-Y')(y-Y'y 

(120) 



Cette expression ne dependant que de E(x,y) et du noyau de Bergmann, elle est totale- 
ment calculable. Cependant, avant de la determiner a la limite conforme, il est utile de 
s'interesser d'abord au cas de W(x,y) dont je n'ai pas encore acheve le calcul. 

5 Remerciements. 

Pour conclure, je tiens a remercier Bertrand Eynard qui, au cours de ces quelques mois 
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